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Abstract 

We present a precise measurement of the standard model CP-violation parameter sin 2(f>i based 
on a sample of 152 x 10 6 BB pairs collected at the T(4S) resonance with the Belle detec- 
tor at the KEKB asymmetric-energy e + e~ collider. One neutral B meson is reconstructed in 
a J/^K° S , %1){2S)K° S , Xc\K° s , r) c K° s , J/ipK* , or J/tpK^ CP-eigenstate decay channel and the 
flavor of the accompanying B meson is identified from its decay products. From the asymme- 
try in the distribution of the time interval between the two B meson decay points, we obtain 
sin 20i = 0.733 ± 0.057(stat) ± 0.028(syst). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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In the standard model (SM), CP violation arises from an irreducible phase in the weak 
interaction quark-mixing matrix [Cabibbo-Kobayashi-Maskawa (CKM) matrix] In par- 
ticular, the SM predicts a CP-violating asymmetry in the time-dependent rates for B° and 
B° decays to a common CP eigenstate fcp, where the transition is dominated by the b — > ccs 
process, with negligible corrections from strong interactions 0|: 

A{t) = T(B°^f CP ) + T(B°^f CP ) = -^ 2 ^HM. (1) 

where T(B°, B° — > fcp) is the rate for B° or B° to fcp at a proper time t after production, 
£/ is the CP eigenvalue of fcp, Am^ is the mass difference between the two B° mass 
eigenstates, and 0i is one of the three interior angles of the CKM unitarity triangle, defined 
as 0i = 7r — &rg(V t lV t d/V^Vai)- Non-zero values for sin20! have been reported by the Belle 
and BaBar collaborations 0, 0, El . 

Belle's latest published measurement of sin 20! is based on a 78 fb -1 data sample (data 
set I) containing 85 x 10 6 BB pairs produced at the T(AS) resonance. In this paper, we 
report an improved measurement incorporating an additional 62 fb _1 (data set II) for a total 
of 140 fb _1 (152 x 10 6 BB pairs). Changes exist in the analysis with respect to our earlier 
result . We apply a new proper-time interval resolution function that reduces systematic 
uncertainties in sin20i and also in Amj and lifetime (t b o, t b +) measurements. We introduce 
6-flavor-dependent wrong-tag fractions to accommodate possible differences between B° and 
B° decays. We also adopt a multi-parameter fit to the flavor-specific control samples to 
obtain the resolution parameters and wrong-tag fractions simultaneously. There are other 
improvements in the estimation of background components, which become possible with 
increased statistics. 

The data were collected with the Belle detector [(J at the KEKB asymmetric collider [?| , 
which collides 8.0 GeV e~ on 3.5 GeV e + at a small (±11 mrad) crossing angle. We use 
events where one of the B mesons decays to fcp at time tcp, and the other decays to a 
self-tagging state /tag' which distinguishes B° from B°, at time ttag- The CP violation 
manifests itself as an asymmetry A(At), where At is the proper time interval between the 
two decays: At = tcp — ttag- At KEKB, the T(45) resonance is produced with a boost of 
/?7 = 0.425 nearly along the z axis defined as anti-parallel to the positron beam direction, 
and At can be determined as At ~ Az/(j3~f)c, where Az is the z distance between the 
fcp and ^ag decay vertices, Az = z C p — ^tag- The average value of Az is approximately 
200 //in. 

The Belle detector Q is a large- solid- angle spectrometer that includes a silicon vertex 
detector (SVD), a central drift chamber (CDC), an array of aerogel threshold Cerenkov coun- 
ters (ACC), time-of-flight (TOF) scintillation counters, and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located outside of the coil is instrumented 
to detect K L mesons and to identify muons (KLM). 

We reconstruct B° decays to the following CP eigenstates [8J]: J/^K^, ip(2S)Kg, XciKg, 
r] c Kg for £f = —1 and J/ipK^ for £j = +1. We also use B° — > J/ipK*° decays where 
K*° — > Kg7f°. Here the final state is a mixture of even and odd CP, depending on the 
relative orbital angular momentum of the J/ip and K*° . We find that the final state is 
primarily £/ = +1; the £/ = -1 fraction is 0.19 ± 0.02(stat) ± 0.03(syst) @- 

The reconstruction and selection criteria for all fcp channels used in the measurement 
are described in detail elsewhere jH. J/ijj and ip(2S) mesons are reconstructed via their 
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decays to (£ = p, e). The ip(2S) is also reconstructed via J/ipir + 7i~, and the xa 

via J/ijyy. The r] c is detected in the K s K~n + , K + K~tt°, and modes. For the J/ipKg 
mode, we use Kg — > 7r + 7r~ and 7r°7r° decays; for other modes we only use — > 7r + 7r~. 
For reconstructed B — > fcp candidates other than J/ipK L , we identify B decays using 
the energy difference AE = E B JRS — ^^ s m and the beam-energy constrained mass M^ c = 

V^^befnn-' 2 — G°B mS ) 2 ' wnere ^beam * s ^ e Deam ener gy i n t ne center-of-mass system (cms) 
of the T(4S) resonance, and i^ ms and p B ms are the cms energy and momentum of the 
reconstructed B candidate, respectively. 

Candidate B° — > J/ipK® decays are selected by requiring ECL and/or KLM hit patterns 
that are consistent with the presence of a shower induced by a K\ meson. The centroid of 
the shower is required to be within a 45° cone centered on the K® direction inferred from 
two-body decay kinematics and the measured four-momentum of the J/ip- 

We perform a multi-parameter fit to flavor-specific control samples to obtain wrong- 
tag fractions and parameters for the resolution function simultaneously. We select B° — > 
D*-£ + u, J/ipK*°(K + n-), D*-n + , D'tt^, D*~p + , and J/tpK s {£ + £-) (for resolution param- 
eters only) for B° decays, and B + — > D°tt + and J j %pK + for B + decays. The total numbers 
of candidates (N ev ) and purities (p) are N e y = 124118 and p = 0.82 for B° decays, and 
iVev = 57305 and p = 0.81 for B + decays. The fit uses free parameters for wrong-tag frac- 
tions (12), for the resolution function (14), for the B + background in B° decays (3), Am<j, 
Tgo and tb+- The total number of parameters in the fit is 32. We add two parameters to 
the resolution function described in ^(j to obtain an improved description of the effect of 
charmed particle decays in the /tag ver t ex - We test the new fit method and parameteri- 
zation with a large number of Monte Carlo (MC) events. A fit to the MC control sample 
yields Am d = (0.488 ± 0.002) ps" 1 , r B o = (1.539 ± 0.003) ps and r B+ = (1.679 ± 0.004) ps 
for the input values of Am d = 0.489 ps -1 , t b o = 1.541 ps and t b + = 1.674 ps, respectively. 
The obtained wrong-tag fractions are also found to be correct. The unbinned maximum- 
likelihood fit to data yields Am d = [0.511 ± 0.005(stat)] ps" 1 , t b o = [1.533 ± 0.008(stat)] 
ps and r B + = [1.634 ± O.Oll(stat)] ps, where the errors are statistical only. The results are 
consistent with the present world average values 

Charged leptons, pions, kaons, and A baryons that are not associated with a reconstructed 
CP eigenstate decay are used to identify the 6-flavor of the accompanying B meson. The 
tagging algorithm is identical to the one used in reference 0. We use two parameters, q 
and r, to represent the tagging information. The first, q, has the discrete value +1 (— 1) 
when the tag-side B meson is likely to be a B° (B°), and the parameter r is an event-by- 
event Monte Carlo-determined flavor-tagging dilution parameter that ranges from r = for 
no flavor discrimination to r = 1 for an unambiguous flavor assignment. It is used only 
to sort data into six intervals of r, according to estimated flavor purity. We determine 
directly from data the average wrong-tag probabilities, wi = (wf + wf)/2 (I = 1,6), and 
differences between B° and B° decays, Awi = — wf , where is the wrong-tag 

probability for the B°(B°) decay in each r interval. The event fractions and wrong-tag 
fractions are summarized in Table |U The total effective tagging efficiency is determined to 
be e g g = J2i=i e «(l — 2wi) 2 = 0.287 ± 0.005, where e\ is the event fraction for each r interval. 
The error includes both statistical and systematic uncertainties. 

The vertex position for the fcp decay is reconstructed using leptons from J/ip decays or 
charged hadrons from r] c decays, and that for /tag ^ s obtained with well reconstructed tracks 
that are not assigned to fcp- Tracks that are consistent with coming from a Kg — > 7r + 7r~ 
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TABLE I: The event fractions q, wrong-tag fractions W[, wrong-tag fraction differences Awi, and 
average effective tagging efficiencies = e/(l — 2wi) 2 for each r interval. The errors include 
both statistical and systematic uncertainties. The event fractions are obtained from the J/ipKo 
simulation. 



1 


r interval 


e« 


Wl 


Awi 


e eff 


1 


0.000 - 


0.250 


0.398 


0.464 ± 0.006 


-0.011 ±0.006 


0.002 ± 0.001 


2 


0.250 - 


0.500 


0.146 


0.331 ± 0.008 


±0.004 ± 0.010 


0.017 ±0.002 


3 


0.500 - 


0.625 


0.104 


0.231 ±0.009 


-0.011 ±0.010 


0.030 ± 0.002 


4 


0.625 - 


0.750 


0.122 


0.163 ± 0.008 


-0.007 ± 0.009 


0.055 ± 0.003 


5 


0.750 - 


0.875 


0.094 


0.109 ±0.007 


±0.016 ± 0.009 


0.057 ± 0.002 


6 


0.875 - 


1.000 


0.136 


0.020 ±0.005 


±0.003 ± 0.006 


0.126 ±0.003 



decay are not used. Each vertex position is required to be consistent with the interaction 
region profile, determined run-by-run, smeared in the r-<p plane to account for the B meson 
decay length. With these requirements, we are able to determine a vertex even with a single 
track; the fraction of single-track vertices is about 10% for zqp and 22% for % a g. The 
proper-time interval resolution function i? s jg(At) is formed by convolving four components: 
the detector resolutions for z C p and ^tag> the sm ^ m ^ ne z tag ver tex position due to sec- 
ondary tracks originating from charmed particle decays, and the kinematic approximation 
that the B mesons are at rest in the cms |l0j]. A small component of broad outliers in the Az 
distribution, caused by mis-reconstruction, is represented by a Gaussian function. We deter- 
mine fourteen resolution parameters from the aforementioned fit to the control samples. We 
find that the average At resolution is ~ 1.43 ps (rms). The width of the outlier component 
is determined to be (39 ±2) ps; the fractions of the outlier components are (2.1 ±0.6) x 10~ 4 
for events with both vertices reconstructed with more than one track, and (3.1 ±0.1) x 10~ 2 
for events with at least one single-track vertex. 

After flavor tagging and vertexing, we find 5417 events in total in the signal region, 
which are used for the sin20i determination. Table HD lists the numbers of candidates, Aev, 
and the estimated signal purity for each fcp mode. Figure [T] shows the M^ c distribution 
after applying mode-dependent requirements on AE for all B° candidates except for B° — > 
J/i/jK^- There are 3085 entries in total in the signal region denned as 5.27 < M^ c < 5.29 
GeV/c 2 . Figure 121 shows the p^ ms distribution for B° —> J/ipK^ candidates. We find 2332 
entries in the 0.20 < p^ ms < 0.45 GeV/c signal region. 

Figure El shows the observed At distributions for the q^f = ±1 and q£f = —1 event 
samples (top), the asymmetry between two samples with < r < 0.5 (middle) and with 
0.5 < r < 1.0 (bottom). The asymmetry in the region 0.5 < r < 1.0, where wrong-tag 
fractions are small as shown in Table U clearly demonstrates large CP violation. 

We determine sin 20i from an unbinned maximum-likelihood fit to the observed At dis- 
tributions. The probability density function (PDF) expected for the signal distribution is 
given by 

e -\At\/r B o 

V s i g {At,q,Wi,Awi,£ f ) = ^ - — [1 - qAwi - q£ f (l - 2wj)sin20i sin(Am d At)], (2) 

where we fix the B° lifetime r^o and mass difference Am^ at their world average values . 
Each PDF is convolved with the appropriate -R s jg.(At) to determine the likelihood value for 
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TABLE II: The numbers of reconstructed B — > fcp candidates after flavor tagging and vertex 
reconstruction, N ev , and the estimated signal purity, p, in the signal region for each fcp mode. J/tp 
mesons are reconstructed in J/tp — > or e + e~ decays. Candidate Kg mesons are reconstructed 

in Kg — > 7r + 7r~ decays unless otherwise written explicitly. 



Mode 


if 


AW 


P 




-l 


1997 


0.976 ± 0.001 


J/ijjK s (Tr Q ir°) 


-l 


288 


0.82 ±0.02 


ip(2S)(£ + e~)K° s 


-l 


145 


0.93 ±0.01 


^{2S){J/^ + ^-)K s 


-l 


163 


0.88 ±0.01 


Xci(J/^)K° s 


-l 


101 


0.92 ±0.01 


r) c (K° s K-ir+)K° s 


-l 


123 


0.72 ±0.03 


7] c (K + K-tt°)K% 


-l 


74 


0.70 ±0.04 


Vc(PP)K% 


-l 


20 


0.91 ±0.02 


All with £ f = -1 


-l 


2911 


0.933 ± 0.002 


J/^K*°(K° s tt°) 


+1(81%) 


174 


0.93 ±0.01 


J/^Kl 


+1 


2332 


0.60 ±0.03 




(GeV/c") 



FIG. 1: The beam-energy constrained mass distribution within the AE signal region for all fcp 
modes other than J/tpK®. The solid curve shows the fit to signal plus background distributions, 
and the dashed curve shows the background contribution. 
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FIG. 2: The p£ ms distribution for B° — > J/i/jK^ candidates with the results of the fit. The dashed 
lines indicate the signal region (0.20 < p^ ms < 0.45 GeV/c). 



each event SIS db function of sin20i: 



Pi = (1 - f \) J [isig^sig( A ^ 9> wi> A ^/, 0)^sig( At " A 

+ (1 - /sig)^bkg( At/ )^bkg( At - At')]rfAt' + / ol P ol (At), 



(3) 



where / s jg is the signal fraction calculated as a function of p^ 1118 for J/ ipK® and of AE and 
for other modes. ^bkg(A£) ^ the PDF for combinatorial background events, which 
is modeled as a sum of exponential and prompt components. It is convolved with a sum 
of two Gaussians, R\^g, which is regarded as a resolution function for the background. To 
account for a small number of events that give large At in both the signal and background, 
we introduce the PDF of the outlier component, P Q i, and its fraction / j. The only free pa- 
rameter in the final fit is sin 20! , which is determined by maximizing the likelihood function 
L = n« Pi, where the product is over all events. The result of the fit is 



sin . 



0.733 ± 0.057(stat) ± 0.028(syst). 



Sources of systematic error include uncertainties in the flavor tagging (0.014), in the 
vertex reconstruction (0.013), in the background fractions for B° J/ipK® (0.012) and 
for other modes (0.007), in the resolution function (0.008), a possible bias in the sin 20! 
fit (0.008), and an effect of interferences [l2| in the / tag final state (0.008). The errors 
introduced by uncertainties in Am d , r B o and in the background At distribution, are less 
than 0.005. 

Several checks on the measurement are performed. Table ITTT1 lists the results obtained by 
applying the same analysis to various subsamples. All values are statistically consistent with 
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-o-q^,=-1 




-2 2 
At(ps) 



FIG. 3: The At distributions for the events with q£f = — 1 (open points) and q£f = +1 (solid 
points) with all modes combined (top), the asymmetry in each At bin between q£f = — 1 and 
q£f = +1 samples with < r < 0.5 (middle), and with 0.5 < r < 1 (bottom). The results of the 
global unbinned maximum-likelihood fit (sin20i = 0.733) are also shown. 



each other. Figure 0] shows the raw asymmetries and the fit results for (cc)K^ (top) and 
J/i)Kl (bottom). A fit to the non-CP eigenstate modes B° -> D*~i + v and J/ipK^iK+Tr"), 
where no asymmetry is expected, yields 0.012 ± 0.013(stat). 

The signal PDF for a neutral B meson decaying into a CP eigenstate [Eq. (J5J] can be 
expressed in a more general form as 



V sig (At,q,w h Awi) 



e -\At\/r B o 

4r s o 



qAwi + q(l — 2wi) S sin(Am^At) + ^4cos(AmrfAt) 



(4) 



where S = 2Xm(A)/(|A| 2 + 1), A = (|A| 2 — 1)/(|A| 2 + 1), and A is a complex parameter 
that depends on both B°-B° mixing and on the amplitudes for B° and B° decay to a CP 
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TABLE III: The numbers of candidate events, 


A ev , and values of sin 2^ 


for various subsamples 


(statistical errors only). 






Sample 


AT 

ivev 


sin 2(p\ 


T Ink Z^O (TT+ir~\ 

J/iPK° s (ir ir°) 
Xcl-^s 


1 QQ7 

ryy i 

288 

OuO 

1 m 


U.O 1 It U.uo 

0.72 ±0.20 
n so -i- n 9n 

i 1^4 _i_ n 40 

J. . (Jt: _l_ W.^i? 




217 


1.32 ± 0.29 


All with £f = — 1 


2911 


73 ± 06 


T/ibK9 

J/i>K* (K° Q -K ) 


2332 
174 


f) 80 ± 1 3 

0.10 ± 0.45 


h ~ — B° (a — +1) 


2717 


72 ± 09 

\J • 1 — _1_ U. Ui7 


f f = B° (q = -1) 


2700 


0.74 ± 0.08 


< r < 0.5 


2985 


0.95 ±0.26 


0.5 < r < 0.75 


1224 


0.68 ±0.11 


0.75 < r < 1 


1208 


0.74 ± 0.07 


data set I (78 fb" 1 ) 


3013 


0.73 ± 0.07 


data set II (62 fb" 1 ) 


2404 


0.74 ± 0.09 


All 


5417 


0.733 ± 0.057 




FIG. 4: The raw asymmetries for (cc)Kg (£j = —1) (top) and J/ipK^ (£j = ±1) (bottom). The 
curves are the results of the global unbinned maximum-likelihood fit. 
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eigenstate. The presence of the cosine term (|A| 7^ 1) would indicate direct CP violation; the 
value for sin 20i reported above is determined with the assumption | A| = 1, as | A| is expected 
to be very close to one in the SM. In order to test this assumption, we also performed a 
fit using the expression above with acp = — £/ImA/|A| and |A| as free parameters, keeping 
everything else the same. We obtain 

|A| = 1.007 ±0.041(stat), 
a CP = 0.733 ± 0.057(stat), 

for all CP modes combined. This result is consistent with the assumption used in our 
analysis. 
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